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ABSTRACT 

Context. The overall evolution of low-mass stars is heavily influenced by the processes occurring in the stellar interior In particular, 
mixing processes in convectively unstable zones and overshooting regions affect the resulting observables and main sequence lifetime. 
Aims. We aim to study the effects of different convective boundary definitions and mixing prescriptions in convective cores of low- 
mass stars, and to discriminate the existence, size, and evolutional^ stage of the central mixed zone by means of asteroseismology. 
Methods. We implemented the Ledoux criterion for convection in our stellar evolution code, together with a time-dependent diffusive 
approach for mixing of elements when semiconvective zones are present. We compared models with masses ranging from 1 Mq to 2 
Mq computed with two different criteria for convective boundary definition and including different mixing prescriptions within and 
beyond the formal limits of the convective regions. Using calculations of adiabatic oscillations frequencies for a large set of models, 
we developed an asteroseismic diagnosis using only / = and Z = 1 modes based on the ratios of small to large separations roi and rio 
defined by Roxburgh & Vorontsov (2003). We analyzed the sensitivity of this seismic tool to the central conditions of the star during 
the main sequence evolution. 

Results. The seismic variables tqi and rio are almost linear in the expected observable frequency range, and we show that their slope 
depends simultaneously on the central hydrogen content, the extent of the convective core, and the amplitude of the sound-speed 
discontinuity at the core boundary. By considering about 25 modes and an accuracy in the frequency determinations as expected from 
the CoRoT and Kepler missions, the technique we propose allows us to detect the presence of a convective core and to discriminate 
the different sizes of the homogeneously mixed central region without the need of a strong a priori for the stellar mass. 

Key words. Stars: evolution ~ Convection - Asteroseismology - Stars: oscillations 
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• 1. Introduction ing processes that has been extensively studied is semiconvec- 

■ tion. During the hydrogen-burning phase the convective core 

O ■ Understanding the physical processes dominating stellar inte- retreats leaving behind a non-uniform chemical profile deter- 

^ riors and correctly characterizing their impact on the evolu- mined by the composition of the core at the moment each 

tion of stars is one of the main goals of stellar astrophysics, layer is detached from it. Depending on the criterion consid- 

> For a long time it has been acknowledged that different mix- ered for defining convective boundaries, these layers could main- 

^ ing processes in stars such as convection, microscopic diffusion tain their chemical composition or undergo some mixing pro- 

and rotational mixing strongly influence the overa ll evolu tion cess, k nown as semiconvective mixin g. Since the pioneering 

of stars (e.g. Zahn 1992; Maeder & Mevnetll2000t iXhoul et al.i work of lSchwarzschild & Hlmildlli^ . several authors have in- 

|1994 iTurcotte et al.,;I998, and references therein), with subse- vestigated the occurrence of se miconvec tive mixing and its ef- 

quent effects in observable quantities such as luminosity, effec- fects on stellar evolution (e.g IStothersl [T970: Stoth ers & Ciiii] 



X 



live temperature and surface abundances. This has a direct im- 1 975 ; Ch iosij 1 978t iLanger et all 1 985TAlthough it was initially 

pact on the determinations of derived parameters relying still to thought to occur only in massive stars, semiconv ection was also 

a large extent on evolutionary models, as in the case of stellar predicted to take place in low-mass stars (e.g. lMitaia^[l972l: 

masses and ages. iFaulkner & CannonI 119731; ISaiol 119751; iGabriel & Noelsf Il977l) 



The treatment of convectively unstable zones and overshoot- but with a larger impact during the helium-burning phase, 
ing regions in evolutionary calculations severely affects the in- During the main sequence evolution, stars with masses larger 

terior structure of the stellar models. This depends on how the than ~1.1 Mq develop a convective core whose extent is deter- 

boundaries of the convective regions are defined and convective mined by their temperature stratification. As the timescale for 

zones themselves treated from the point of view of mixing of mixing of elements in this convective region is much shorter 

elements. The resulting luminosity, temperature, and main se- than the nuclear timescale, the core is believed to be homoge- 

quence lifetime of the star are largely affected by these pro- neously mixed and a discontinuity in density appears at the edge 

cesses. In the case of massive star evolution, one of the mix- of the fully mixed core. This discontinuity is produced either by 
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convective core expansion due to the increasing importanc e of 
the CNO cycle over the pp chain ("Mitalas HlQT^ISaiolfTOTSl) . or 
by the retreating convective core leaving behind a non-uniform 
chemical profile (iFaulkner & CannonI 19731) . Both cases produce 
higher opacities outside the convective core and allow a semi- 
convective region to develop. If the resulting density barrier is 
sustained throughout the main sequence evolution the homo- 
geneously mixed central region will be either restricted in its 
growth or not allowed to develop at all, drastically changing 
the behavior of the evolutionary tracks especially close to the 
main sequence turn-off. This has important implications for the 
use of the Color-Magnitude Diagram (CMD) morphology and 
the existence of a hook-like feature at the turn-off to estimate 
properties of stars at the end of the main seq uence and the ag e 
of a given stellar population (see for instance lMaededll974allbl: 
,Maeder & Meynet 1991). 

While it is not possible to obtain direct information about 
stellar interiors, helioseismology has provided the best example 
of indirect observations from the interior of a star by piercing the 
outer layers of the Sun using the oscillations observed at its sur- 
face. These data have served to extract the adiabatic sound speed 
and density profiles of the Sun, determine the location of the so- 
lar helium second ionization zone and the base of the convective 
enve lope, as well as to constrain the solar surfa ce abundances 
(see IChristensen-D alsgaardl f2002l: iBasu & Antial 1200 8. and ref- 
erences therein). With this principle in mind, efforts are being ad- 
dressed to extend the techniques applied to the Sun to other stars 
by means of asteroseismology, i n the hght of the d ata currently 
being obtained by the CoRoT jBagUn et al.ll2006h and Kepler 
(IBorucki et al.ll200 9) missions. 

We address in this paper the issue of the determination of 
convective boundaries and the treatment of zones which present 
a gradient in the molecular weight by the inclusion of the Ledoux 
criterion for convective instability, and a diffusive approach for 
semiconvective mixing in our stellar evolution code. We explore 
the possibility of determining the presence of a convective core 
and discriminating different sizes of homogeneously mixed cen- 
tral regions and evolutionary stages by means of asteroseismol- 
ogy. We do this by exploring the effects that different types of 
convective mixing processes have on asteroseismic frequency 
combinations, in agreement with the expected quality of data 
being obtained by the aforementioned space missions. 

The paper is organized as follows. In Sect. |2]we describe 
the differences in the definition of the convective boundaries for 
each stability criterion, and the prescription for mixing in con- 
vective, semiconvective and overshooting regions. Sect. |3] fo- 
cuses in the impact of these boundary definitions and mixing 
processes during the main sequence evolution. In Sect.|4]we dis- 
cuss the method of isolating the influence of the stellar core by 
means of different asteroseismic tools. We assess the possibility 
of disentangling between the different mixing prescriptions us- 
ing the expected observable range of frequencies. A discussion 
on the open issues and detection possibilities is given in Sect. |5] 
We conclude and state final remarks in Sect. |6] 



2. Stellar Modeling 

We used the Garchins 



Stellar Evolution Code (GARSTEC, 

IWeiss & Schlattll l2008l) for our calculations. The interested 
reader can find in that reference a more detailed description 
of the numerics and input physics included in the code. For 
our model calcul ations we used the 2005 version of the OPAL 
equation of state (iRogers et al.lll996tlRogers & Navfonovl2002l) 
complemented with the MHD equation of state for low tem- 



peratures (iHummer & Mihalas 1 19881) . low temperature opaci- 
ties from iFerguson et al.' (2005^ and OPAL opaciti es for high 
temperatures (llgle sias & Rogers 1996), the iGreves se & Sauva| 
(Il998h solar mixture, and the NACRE compilation for ther- 
monuclear reaction rates (lAngulo et al.,il999) . Convective zones 
are treated with the mixing-length theory (MLT), no matter 
which criterion is used to define their boundaries. Within these 
zones, the chemical composition is modified either instanta- 
neously or by a diffusive process using the convective ve- 
locity es timated from t he ML T as described, for instance, in 
Kippenh ahn & WeigertI ([T990h . 

For the frequency computations of the specific models an- 
alyzed, we have used the Aarhus adiab atic oscillation package 
rADIPLS. IChristensen-Dalsgaardl2008l) . 

2. 1 . Boundaries of convective and semiconvective zones 

The convective zones within a stellar model are defined 
using a stability criterion. The convectively unstable re- 
gions ar e usually determined by the Schwarzschild criterion 
( Schwa rzschild & H arm 1958) by comparing the radiative and 
adiabatic temperature gradients, denoted V^d and V^d respec- 
tively. A layer is convective when 



V,,d < V„ 



(1) 



However, gradients of molecular weight p can have a stabilizing 
- or desta bilizing - effec t on the convection process. The Ledoux 
criterion (lLedouxlll947l) takes into account this effect. Within 
this prescription, a layer is convective when 



Vl < V„d, 

with Vl, the Ledoux temperature gradient, defined as 
Vl = Vad + 



where 



(dlnp),/^ (sinrjp/^'' ( 



dlnp 
dlnP 



(2) 



(3) 



(4) 



For an equation of state appropriate to a mixture of an ideal gas 
and black body radiation, Eq. |3]reduces to: 



Vad + 



/3 



(5) 



where 13 is the ratio of gas pressure to total pressure 
(iKippenhahn & Weigertlll990l) . 

The inclusion of changes in the molecular weight in the def- 
inition of the relevant temperature gradients leads to the ap- 
pearance of zones whose energy transport process will depend 
on the convective criterion considered. There are regions which 
would be considered convective if the Schwarzschild criterion 
were used, but are in turn radiative if the Ledoux criterion were 
applied. These zones are called semiconvective zones, and are 
defined as the regions where 



Vad < Vrad < Vl . 



(6) 



As an example, a model of a 1.5 Mq star displaying a semi- 
convective zone during its main sequence evolution is depicted 
in Fig.[T] The hydrogen profile in the stellar interior is constant 
inside the homogeneously mixed convective core, and presents 
a sharp variation at the position where the edge of the core is 
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Fig. 1. Internal structure of a 1.5 Mq stellar model in the main 
sequence. The hydrogen abundance profile X (black solid line), 
the adiabatic gradient Vad (red dotted line), the radiative gradient 
Vrad (green dash-dotted line), and the Ledoux gradient Vl (blue 
dashed line) are plotted as a function of the mass fraction. From 
left to right, tones of gray fill regions of different energy transport 
processes: convective zone, semiconvective zone, and radiative 
zone. 

located. This feature produces both a step in the radiative gra- 
dient Vrad due to the change in the opacities, and the departure 
of the Ledoux gradient Vl from the adiabatic gradient Vad as a 
result of the change in the molecular weight. Therefore, a signif- 
icant semiconvective layer appears at the boundary of the con- 
vective core. We notice that, in this example, the semiconvective 
zone is 60% as massive as the convective region itself. If the 
Schwarzschild criterion were used instead, the semiconvective 
layer would be part of the convective core. 

2.2. Prescription for semiconvective layers 

Throughout the years, several approaches have been proposed 
to deal with semiconv ective r egions when they appear in stel- 
lar models ([St evenso n [p79l iLanger et alJ [T98l ISpruid [T99l 
iGrossman & Taamlll99q) . According to the used prescription, 
the semiconvective zone is considered to be mixed more o r less 
efficiently (see the comparison made by Merrvfieldll 1 9951) . The 
amount of mixing in semiconvective layers is crucial: an efficient 
mixing process can reduce the molecular weight gradient suffi- 
ciently, and then trigger convective instability. For the present 
work, we ha ve imp lemented in GARSTEC the prescription pro- 
posed bv lLanger et al.l (Il983 . 1985), based on the description of 
vibrational instability made by .Kato, (.1966) . With this method, 
the true temperature gradient V in a semiconvective region is 
calculated from the relation 



imd °'"2V(Vl-V) 



(V - Vad) - 



/^(8-3/j) , 
32-24^-/?2 



(7) 



where Lsc and L^d are the semiconvective and radiative lumi- 
nosities, respectively. 

The mixing is treated as a time-dependent diffusive process, 
with the diffusion coefficient calculated as 



^sc - Q's( 



v-v 



where is the radiative conductivity, Cp the specific heat at 
constant pressure, and p the density. Both the diffusion co- 
efficient and the temperature gradient in this prescription de- 
pend on an efficiency parameter of semiconvection (asc)- The 
meaning of this parameter is further discussed in Sect. 13.2] 
We stress the point that this diffusive and time-dependent ap- 
proach significantly differs from the ones used in previous ef- 
forts to treat semiconvection in low-mass stars. In previous 
studies, the mixing in the semiconvective layer is either per- 
formed by adjusting the composition until convective neutral- 
ity according to the Schwarzschild criterio n is again reached 
^owe & Mitalas 1982; Miglio et al. 2008), or approximating 
the results obtained by Spruit ( 1 992.) with a two-step function 
dPopielski & Dziembowskill2005h~ 

2.3. Oversliooting 

In some of our models we have considered mixing of chemical 
elements beyond the formal convective boundaries. We assume 
that this process induces mixing but does not modify the thermal 
structure of the layer s (overshooting, if we use the terminology 
given bv IZahnlll991h . It is implemented in our code as a dif- 
fusive process consisting of an exponential de chne of the con - 
vective velocities within the radiative zone (Frevtag et al.lll996h . 
The diffusion constant is given by 



Dov (z) = Do exp 



-2z 



(9) 



ad 



6 Cp p Vl - V ' 



(8) 



where ^ corresponds to an efficiency parameter. Hp is the pres- 
sure scale height, z is the distance from the convective border, 
and the diffu sion constant Dq is derived from MLT-convective 
velocities ( Ki ppenhahn & Weigertlll990l) . We have used an over- 
shooting efficiency of f = 0.016, which is in t he range expected 
for main sequence stars (iHerwig et al.l 1 19971) and corresponds 
to the value obtained by calibrating the parameter ^ with open 
clusters. The size of the overshooting region is further limited in 
the case of small convective cores. This is done in GARSTEC 
using a geometrical cutoff factor, allowing the overshooting re- 
gion to e xtend only to a fraction of the convective zone (see 
IWeiss & Schlattl.2008; .Magic et aL.2010. for details). 

3. Impact on the Stellar Models 

Several models were computed using the different mixing pre- 
scriptions and definition of the convective boundaries. Within a 
convective zone the mixing is performed diffusively using the 
MLT while in a semiconvective region the mixing is carried out 
as explained in Sect. 12.21 Calculations were made for both the 
Schwarzschild and the Ledoux criterion for the definition of the 
convective zones, with and without including extra mixing due 
to overshooting and semiconvection. All the models considered 
here are computed for solar metallicity. We explore the effects 
of these processes for masses ranging from 1 to 2 M0 starting 
at the pre-Zero Age main sequence and evolved until hydrogen 
exhaustion in the core. 



3.1. Convective boundary definition criteria 

Applying either the Schwarzschild or the Ledoux criterion in- 
fluences the appearance and the size of convective cores as a 
function of stellar mass. To study these effects, we computed 
models with no convective overshoot and no mixing throughout 
the semiconvective regions when these were present (ajc - 0). 
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In Fig.|2]we present the evolution during the main sequence 
of convective and semiconvective regions in the interiors of three 
representative models: 

- a model where the inclusion of the Ledoux criterion inhibits 
the growth of the convective core in the late phase of the 
main sequence evolution (1.2 Mq model, Fig.|2]fop panel); 

- a model with a convective core increasing in size during the 
hydrogen burning phase whose growth is constrained if the 
Ledoux criterion is applied (1.5 Mq model. Fig. ^middle 
panel); 

- a model where the convective core recedes during the main 
sequence leaving behind a chemical discontinuity (2.0 Mq 
model, Fig.^bottom panel). 
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Fig. 2. Evolution of the convective core size in mass coordinate 
as a function of central hydrogen content during the main se- 
quence phase for the two considered convective boundary defini- 
tion criteria: Ledoux (solid blue lines) and Schwarzschild (solid 
red lines). Blue dotted lines depict semiconvective zones when 
the Ledoux criterion is applied (by definition there are no semi- 
convective zones when the Schwarzschild criterion is applied). 
The panels represent different masses: 1.2 Mq (top), 1.5 M0 
(middle) and 2.0 M0 (bottom). 



There are clear differences in the convective core behavior for 
the two considered criteria. In the case of growing convective 
cores (Fig. |2] top and middle panels) the molecular weight dis- 
continuity at the edge of the convective zone produces a steep 
increase in the opacity, which translates into a sharp discontinu- 
ity in the radiative gradient (as shown in Fig. [TJ. 

If the Schwarzschild criterion is applied, the profile of the 
convective core presents wrinkles throughout the evolution, sig- 
nature of a convectiv e process not properly taken into account 
(iLebreton et al.ll2008h . The exact position of the convective core 
edge will depend largely in the numerics of the evolutionary 
code and the way the mesh points are place d where the discon - 
tinuity in the molecular weight is located (iMigho et al.ll2008h . 
Slight variations in the radiative gradient are induced by, for in- 
stance, the interpolation in the opacity tables, the allowed time- 
step and the accuracy of the calculations, resulting in a small 
region right outside the core edge becoming convective and sup- 
plying fresh hydrogen to the outer layer of the core. The core 



expands and produces the crumple profiles (backwards loops) 
observed in the top and middle panels of Fig. |2] in a phenom- 
ena similar to the so- called breathin g pulses in horizontal branch 
stars (see for example [Catelanl2007l for a detailed discussion and 
further references on the topic). 

If the Ledoux criterion is applied instead, the convective 
core profile remains flat and continuous throughout the main se- 
quence evolution, as the molecular weight discontinuity at its 
edge inhibits its growth. The semiconvective zone appearing out- 
side the core is not mixed through and thus becomes larger dur- 
ing the hydrogen burning phase. At the end of the main sequence 
evolution, the convective core recedes leaving a radiative layer 
between it and the semiconvective zone, causing a change in the 
temperature gradients and making the semiconvective region to 
ultimately disappear. 

In a model with a receding core the convective central re- 
gion reaches its largest size during the initial phase of hydrogen- 
burning, the Schwarzschild criterion allowing it to grow to a 
slightly larger extent than the Ledoux criterion (Fig. |2] bottom 
panel). The chemical discontinuity left behind by the retreat- 
ing core leaves an imprint in the molecular weight which per- 
mits a semiconvective zone to develop, but the radiative gradient 
decreases following the shrinking of the core transforming the 
zone into a radiative one. Since the effect is almost negligible for 
this case, we will focus for the rest of the paper on the models 
where the inclusion of the Ledoux criterion plays an important 
role (masses below ~1.7 M0 with growing convective cores). 

3.2. Main sequence evolution 

In the previous section we have shown the impact that differ- 
ent mixing prescriptions have on the size of the homogeneously 
mixed central region during the main sequence evolution. The 
size (and existence) of the convective core will determine the 
appearance (or lack) of a hook-like feature in the CMD at the 
end of the hydrogen-burning phase. This feature at the turn-off 
is used in the determination of the age of stellar populations by 
fitting isochrones reproducing the shape of the hook at a cer- 
tain metallicity, thus setting constraints on the stellar mass at the 
end of the main sequence for a given set of input physi cs (as 
an example, see IVandenBerg et al.ll2007l: iMagic et alj|201 0, for 
the specific case of M67). As discussed in Sect. 13.11 using the 
Ledoux criterion can inhibit the growth of the convective core 
for the 1.2 Mq model and therefore the appearance of the hook 
feature in the evolutionary track. This is the case when no mix- 
ing is performed within the developing semiconvective region 
(ffsc - 0) and the molecular weight barrier caused by hydrogen 
burning does not allow the core to grow. 

The amount of semiconvective mixing is controlled by the 
efficiency parameter a^c and has a value restricted to a^c < 1, 
but has been used usually in the range 0.001 < ffsc < 0.1 
(iLanger etal.l [l985: Langer 1991; Merryfield 1995). Higher val- 
ues of a^c mean faster mixing velocities, achieving instantaneous 
mixing when a^c 00. In Fig. [3] we have plotted evolution- 
ary tracks for 1.2 and 1.5 M0 cases calculated with different 
mixing prescriptions and semiconvective efficiency parameters. 
By increasing the value of ffsc the convective core develops to a 
larger extent and the evolutionary tracks resemble the one calcu- 
lated with the Schwarzschild criterion. For the 1.2 Mq case the 
track computed with a^c - 0.001 does not increase the size of 
the convective core with respect to the track calculated with no 
mixing in the semiconvective zone, thus their tracks overlap in 
the Hertzprung-Russel Diagram (HRD). One interesting aspect 
is that a position in the HRD can be shared by stars with the 



V. Silva Aguirre, et al.: Probing mixing in stellar cores using asteroseismology 



5 



same mass but different evolutionary stages, internal structures 
and stellar ages. 
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Fig. 3. Hertzprung-Russell Diagram for 1.2 M0 {left panel) 
and 1.5 Mq {right panel) models. Computations done with the 
Ledoux criterion and no mixing in the semiconvective regions 
are plotted with dashed blue lines, Schwarzschild criterion with 
solid red lines, and dotted black lines correspond to tracks with 
Schwarzschild criterion and overshooting. Two values for the 
semiconvective efficiency ffjc are also plotted in thin green dash- 
dotted line (ffsc = 0.001; not shown for the 1.2 Mq as it over- 
laps with the dashed track) and thick light blue dash-dotted line 
(csc = 0.01). The tracks span the evolution from the Zero-Age 
main sequence until hydrogen exhaustion in the core. 

Changing the prescription for convection has an impact on 
the stellar age at which the end of the main sequence is reached. 
In the diffusive prescription we use for semiconvective zones, a^^ 
modifies the value of the diffusion coefficient and also the true 
temperature gradient in the semiconvective zone. The growth 
of the convective core is not enhanced but restricted due to the 
molecular weight gradient and no extra fuel supply is added to 
the core, which translates into shorter main sequence lifetimes. 
In Fig. ID we present the hydrogen contents at the center of the 
models as a function of age, where the differences are clearly 
visible for the considered mixing prescriptions and can amount 
to 30% of the main sequence lifetimes. 

3.3. Effects of mixing in semiconvective layers 

In Fig.|5]we present hydrogen profiles near the center for 1.5 Mq 
models calculated with the Ledoux criterion and different semi- 
convective mixing efficiency values, with the Schwarzschild cri- 
terion and with overshooting. The sizes of the homogeneously 
mixed zones are different as well as the shape of the chemical 
profile at the edge of the convective core. It is worth noticing 
that for the 1 .5 Mq case considered in this paper, a semiconvec- 
tive coefficient of o-sc ~ 0.01 is already efficient enough to very 
closely reproduce the results obtained when the Schwarzschild 
criterion is applied (see the evolutionary tracks in Fig. |3). The 
same is true for c^sc ~ 0.1 in the 1.2 Mq model. This suggests 
that the value of the semiconvective efficiency required to repro- 
duce the results obtained with the Schwarzschild criterion de- 
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Fig. 4. Central hydrogen content of the selected models dur- 
ing their main sequence evolution: 1 .2 Mq (blue) and 1 .5 Mq 
(red). Solid lines: Schwarzschild criterion without overshooting. 
Dashed lines: Ledoux criterion with no semiconvective mixing 
(ffsc = 0). Dash-dotted lines: Ledoux criterion and a semicon- 
vective efficiency of ffjc = 0.01. Dotted lines: models with con- 
vective overshooting. 

creases with mass. In relative terms, the size of the semiconvec- 
tive zone compared to the total convective core size is smaller at 
higher masses, leading to a shorter mixing timescale. 

Convective overshooting has for long been thought of as a 
natural process capable of smoothing out molecular w ei ght gra- 
dients close to the conv ective core (e . g. iMaeded 1 974al lbt iLaneeJ 
ll99lHNoels et al.l20Tol) . We computed models using the Ledoux 
criterion and the overshooting prescription described in Sect. 12. 31 
and realized that this is also the case for low-mass stars using 
our calibrated value for the overshooting efficiency. One should 
keep in mind, though, that the calibratio n itself was based on cer- 
tain assumptions about convection (e.g iPietrinfemi et al.ll2004t 
VandenBerg et al...2006.) . 

4. Asteroseismic diagnoses 

We have highlighted in Sect. |3] the importance of the different 
mixing prescriptions considered in this paper on structure and 
evolution of low-mass stars. In practical terms, these convective 
processes translate into different convective core sizes, ages, and 
evolutionary stages for compatible positions in the HRD. The 
quest of detecting the variations in stellar interiors produced by 
differences in the mixing and energy transport processes can be 
fulfilled by a technique capable of piercing the outer layers of 
stars and which is sensitive to density discontinuities. In this re- 
spect, asteroseismology is thought to be the most powerful tool 
to study the inner proper ties of stars through the observation 
of stellar oscillations (e.g 'Cunha et al."2007l: lAerts et al.ll2008l: 
Christensen-Dalsgaard & Houdek 2010, and references therein). 

Many efforts have been conducted to discover the physics 
of stellar interiors by understanding the influence that dif- 
ferent physical phenomena have on the observed frequency 
spectra. In particular, it is known that a discontinuity in the 
chemical profile of a star produces a sharp variation in the 
adiabatic sound speed, which i n turn introduces an oscillatory 
component in the frequencies (IVorontsovl 1 1 9881: lGoughl[l990l: 
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Fig. 5. Hydrogen profiles for similar central hydrogen con- 
tents of the 1.5 M0 models and different mixing prescriptions: 
Ledoux criterion and no mixing in the semiconvective zones 
(dashed blue line), a^^ - 0.001 (thin dash-dotted green line), 
- 0.01 (thick dash-dotted light blue line), Schwarzschild cri- 
terion (solid red line) and one model including convective over- 
shooting (dotted black line). 



iProvost et al.l Il993l) . This fact has already been applied to 
determine the position of the base of the conve ctive envelope 
in the Sun (Christensen-Dalsgaard etjl. Il99l|), t o constrain 
the amount of overshoot bel o w it ^ onteiro et al.l 1994t 
Roxburgh & Vorontsovl Il994al: IChristensen-Dalsgaard et all 



1995h . and to constr ain the properties of the solar helium second 
ionization z one (Christensen-Dalsgaard & Perez Hernandez 
119911; lAntia'&Basu 1994; Monteiro & Thompson 2005). The 
same technique has been proposed to determine the position 
of the base of the convective envelope and estimate the helium 
abundance for stars other th an the Sun (Monteiro et al. 2000; 
iMonteiro & Thompsonlll998h . 

Following this principle, convective cores in low-mass stars 
have been the subject of several studies by means of asteroseis- 
mology aiming to determine the evolutionary state and size of 
the mixed central region in stars (e.g Audard & Provost 1994[^ 
Audard et al. 1995; Mazumdai" & Antiall2001l: iMaz umdar et al. 
2006; .Cunha & Metcalfe .2007,) . Most of these investigations 



have focused their attention on constraining the size and chemi- 
cal composition of the convective core, while only a few of them 
have acknowledged the impact of its boundary definition and ex- 
tra mixing processes such as semiconvection and overshooting 
in the oscilla t ion frequenci es (*Popielski & Dziembowskil l2005t 
lGodartll2007t iMigho et all 12008 ;"Lebreton et al. 2009). In this 
section we aim at finding a seismic tool that isolates the stellar 
core and which is sensitive to its size and central hydrogen con- 
tent, allowing to disentangle between the different mixing pre- 
scriptions. 

4.1. Seismic variables suited to isoiate tiie core 

In order to study the inner structure of a star using asteroseismol- 
ogy, first we need to find an appropriate seismic variable which 
allows us to probe the desired region of the star and extract the 
required information. Different combinations of low-degree p- 
modes have been suggested as suitable probes of the physical 



characteristics of a star (e.g. IChristensen-Dalsgaardlll984l) . the 
most commonly used being the so-called large and small fre- 
quency separations defined as: 



A/(n) = v„j - Vn-\j 

dlj+2(n) = V„j - Vn-\,l+2 : 



(10) 
(11) 



where v„j is the mode frequency of angular degree I and ra- 
dial order n. These combinations are affected by the outer layers 
of the star where turbulence in the near-surface is almost ne ver 
taken into account into star modeling (e.g. Ball ot et al.ll2004l) . 

In our case, we need to isolate from surface contamination 
the signal arising from the interior of the star in order to properly 
quantify the effects in the frequency spectra of the presence of a 
convective core. Roxburgh & Voront sov (2003) proposed to use 
the smooth 5 points small frequency separations and the ratio of 
small to large separations, and showed that these quantities are 
mainly determined by the inner structure of the star. They are 
constructed as: 



dmin) = ^(v„-i,o - 4v„_i,i + 6v„,o - 4v„,i -i- v„+i,o) 
dw(_n) = --^(v,,-!,! - 4v„,o + 6v„,i - 4v„+i,o + v„+i_i) 



?"oi(n) 



do\{n) 
Ai(«)' 



no 



d\dn) 
Ao(n + 1) 



(12) 



(13) 



(14) 



The small frequency separations have already been used to iden- 
tify the location of the convective envelope and the helium sec- 
ond ionization zone in the Sun by applying them to observational 
data (Roxburgh 2009), while it has been shown that the ratios 
fairl y cancel out the influence in the frequencies of the outer lay- 
ers (lRoxburghll200i l2005t lOtf Floranes et al 1 120051) . As men- 
tioned before, sharp variations on the adiabatic sound speed pro- 
duce an oscillatory signal in the p-mode spectra whose period is 
related to the position of such a discontinuity. In particular, the 
period of this oscillation relates to the travel time of the wave 
through the acoustic cavity, and therefor e to the radial coordi- 
nate where the discontinuity i s loca ted dMonteiro et al.l [T994I : 
Roxburg h & Vorontsovl[T9945l2001b . This travel time is repre- 
sented by the acoustic radius, given by: 



r dr 

Jo Cs 



(15) 



where is the adiabatic sound speed. There is an alternative rep- 
resentation of the acoustic radius called the acoustic depth (r), 
which measures the travel time from the surface towards the in- 
terior Therefore, if the total acoustic radius is given by - t(R), 
then clearly t = - f . If the location of the density discontinuity 
is given by, say, ri in radial coordinates, and that same position 
is represented in acoustic radius and acoustic depth by f i and ri 
respectively, the periods of the oscillation in frequency space in- 
duced by this sharp variation are l/(2fi) and l/(2ri), according 
to the considered variable. That being the case, if either the ra- 
tios or the small separations are affected only by the position of 
the convective core, we could extract this information from the 
frequency data. 

4.2. Sensitivity of and roi to tiie core 

If the differences and ratios defined in Eqs. [12] [13] and [14] 
are indeed sensitive to the deep stellar interior, they should be 



V. Silva Aguirre, et al.: Probing mixing in stellar cores using asteroseismology 



7 




300 1900 3500 5100 100 200 300 6150 6000 5850 5700 
Frequency (/^Hz) Acoustic Radius (s) Tgff (K) 



Fig. 6. 1.1 Mq models with a convective core (green) and without one (red). Each row depicts models sharing the same large 
frequency separation. Left Panels frequency ratios as defined in Eq. [14] Horizontal lines show a typical observable frequency range 
for each mixing prescription of radial orders n = 15-27. Central Panels: adiabatic sound speed as a function of the acoustic radius. 
Right Panels: evolutionary tracks for the two cases considered, the position of the studied models being represented by a filled 
circle. Values of the selected large frequency separations are given (in //Hz) 



strongly influenced by the presence of a convective core or 
a steep gradient in the adiaba tic sound speed near the center 
(iRoxburgh & Vorontsovll2007l) . To test this, we produced two 
sets of evolutionary tracks for a 1 . 1 Mq star, one using the nor- 
mal Schwarzschild criterion and the other one including over- 
shooting without the geometrical restriction described in Sect. 
12.31 This way, we produced a 1 . 1 Mq model with a growing con- 
vective core during the hydrogen burning phase. In order to com- 
pare the effects of the presence of the core, we plotted the evo- 
lution during the main sequence of these models in Fig. |6] The 
left panels present the ratios as a function of frequency, while the 
center panels show the adiabatic sound speed near the center for 
the selected models, whose position in the observational plane is 
depicted by a filled circle in the right panels of the figure. 

The frequency ratios present a common component in the 
low frequency domain (see left panels), which therefore should 
be determined by the global seismic properties of the star: the 
selected models have the same large frequency separation. At 
higher frequencies, the ratios deviate from each other revealing 
the structural differences between the models. This part of the 
frequency domain presents two contributions: one given by the 
presence of a convective core and another one dominated by the 



evolutionary state of the model (equivalent to the central hydro- 
gen content). The oscillatory component observed in the ratio of 
the model calculated with overshooting is related to the presence 
of a convective core. 

A straight line can be drawn between the point where the 
ratios start deviating from each other and the frequency value 
where the oscillatory component of the model with the convec- 
tive core shows its first minimum. We will refer to this part 
of the frequency domain as the linear range from here on. 
Independently of the presence of a convective core or not, we 
observe an increase of the slope - strictly speaking of its abso- 
lute value - with evolution, suggesting that it traces the central 
hydrogen content of the star (and hence stellar age). Thus, the 
sole presence of a negative s lope does not ensure the existence 
of a convective core (see also lBrandao et al.ll2010l) . 

However, in the last row of panels in Fig.|6]the models plot- 
ted have very different central hydrogen abundances (Xc ~ 0.14 
versus Xc ~ 0.4 for the model without and with convective core, 
respectively). In the model without convective core, the increase 
of the slope is mainly due to the growth of the density gradient 
in the center, built during hydrogen burning. Nevertheless, the 
slope is sensibly higher for the model with a convective core. 
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even if this star is less evolved. This difference is a clue to the 
presence or absence of a convective core and will be discussed 
in Sect. 144] 

The global behavior of the small frequency separations is the 
same as the one observed in the frequency ratios. Thus, we will 
focus the analysis in the ratios keeping in mind that the same 
conclusions are applicable to the case of the separations. 

4.3. Influence of mixing prescriptions on rm and roi 

When considering only models which do have a convective 
core, we can study the sensitivity of the ratios to the position 
and size of the convective core, and to the size of the den- 
sity discontinuity. We explore the behavior of 1.3 M© mod- 
els in the three most extreme cases in terms of the size of the 
mixed region: Ledoux criterion with no extra mixing in the 
semiconvective region, the normal Schwarzschild criterion, and 
one evolutionary sequence including overshooting. An error of 
0.2 juHz in the frequency determinations w as assumed, as ex- 
pected from the CoRoT and Kepler missions (iBaglin et al.l2006t 
[Christensen-Dalsgaard et al. 2007). 

In Fig.|7]we present a similar diagram as in Fig. |6]showing 
the changes of this seismic diagnose as a function of the central 
hydrogen content, mapping therefore the main sequence evolu- 
tion. When considering one mixing prescription (color) we no- 
tice that the absolute value of the slope of the frequency ratios 
increases with evolution. The considered models have the same 
central hydrogen content, so when observing this behavior of the 
ratios nothing can be said about the exact evolutionary stage of 
the star However, the size of the discontinuity in the adiabatic 
sound speed is indeed different, as well as the size of the cores, 
the largest being the one with the most negative slope. 

One is tempted to interpret this as a direct relation between 
the size of the discontinuity and the absolute slope of the ratios, 
but two details must be taken into account: for the case of the 
Ledoux criterion, the change in the slope during the main se- 
quence evolution is much smaller than for the other cases, and 
the acoustic radius of the fully mixed region remains almost con- 
stant (while for the other prescriptions it increases). It is clear 
that the value of the absolute slope has contributions coming 
from the position of the convective core and the size of the sound 
speed discontinuity. 

Bearing this in mind, instead of using the central hydrogen 
content as a proxy for the evolution of the ratios we considered 
models with the same large frequency separation, as it is the pa- 
rameter first and most easily obtained from seismic observations. 
Figure [8] shows this case for the 1.3 Mq models with large sep- 
arations equal to Av =110, 100, 95 and 88 yuHz. Due to the dif- 
ferences in the mixing prescriptions, the selected models do not 
have the same central hydrogen content but almost the same age. 

The behavior at the low frequency range (v < 1200 yuHz) 
is the same for the three mixing prescriptions for a given large 
frequency separation value. As already mentioned in Sect. 14.2! 
this is an indication of this part of the frequency spectrum being 
dominated by the global seismic properties of the star, while the 
frequency values above ~ 1200 yuHz seem to be sensitive to the 
differences in the interior (the linear range). 

Within this frequency region where the information of the 
core is contained, we can compare in Fig. |8]a case opposite to 
that in Fig.|7l the last row presents models with the same size of 
discontinuity of the adiabatic sound speed (~ 0.27 x 10^ cm/s). 
The hydrogen content of the Ledoux (blue) model is only ~ 2% 
while for the the Schwarzschild criterion and the model with 
overshooting it is ~ 28% and ~ 32%, respectively. The model 



calculated with the Ledoux criterion again has the smallest ab- 
solute slope value and size of the convective core. 

We stress the fact that the absolute value of the slope cannot 
be attributed only to the size of the adiabatic sound speed discon- 
tinuity. In fact, we can distinguish three effects directly linked to 
the global behavior of the ratios: 

- a decrease in the central hydrogen content produces an in- 
crease in the slope in the linear range, which is visible no 
matter if a convective core is present or not; 

- when a convective core is present, an increase of the size 
of the discontinuity in the adiabatic sound speed at the core 
boundary increases the amplitude of the oscillatory signa- 
ture, making the slope of the ratios steeper; 

- the extent of the convective core also affects the slope: a 
larger convective core reduces the period of the oscillatory 
component in the ratios, increasing in turn the absolute value 
of the slope. 

It is not straightforward to determine how much influence on 
the slope each of these effects have. A differential comparison 
must be done in order to assess where the major contributions 
are coming from. Nevertheless, now that the key processes af- 
fecting the frequency ratios have been identified, we can use this 
information to disentangle between models with different central 
mixed zones sizes and evolutionary stages. 

4.4. Breaking the degeneracy: splitting mixing prescriptions 
and finding convective cores 

In the examples above we have shown that the frequency ratios 
are sensitive to the central hydrogen content, the existence and 
extent of a convective core, and the amplitude of the discontinu- 
ity in the adiabatic sound speed. Thus the importance of being 
able to disentangle between the differe nt mix ing prescription s. 
An attempt to do so was presented bv Popielski & Dziembowskil 
02005) . where they characterized the size of the mixed zone in 
stellar models using the small frequency separation. The authors 
studied differences between two types of semiconvective mixing 
and two values of convective overshooting by averaging over a 
range of radial orders the small frequency separation and using 
its slope. We will apply a similar approach to the frequency ra- 
tios using the radial orders n= 15-27, a typical range which is 
expected to be detected by the CoRoT and Kepler missions (see 
horizontal lines in the left panels of Figs.|6]|7]and[8]l. 

In Fig. |9] we present a slope versus mean diagram in the ob- 
servable range of the frequency ratios, for the complete main se- 
quence evolution of different masses and mixing prescriptions. 
For each mass value and mixing prescription we have also se- 
lected one model with a large frequency separation value of 96 
juHz, their properties presented in Table [T] The position of the 
zero-age main sequence model for a given mass value is the same 
independent of the applied mixing prescription, as the interiors 
are until that point the same. Once the models start evolving, the 
different behaviors in their frequency ratios make them clearly 
distinguishable. The tracks reach a minimum value in slope near 
the exhaustion of the hydrogen in the center (when the convec- 
tive core starts receding if present). 

For the 1 . 1 Mq tracks the two cases present very different 
values of slopes throughout the evolution, allowing us to disen- 
tangle between a star with a convective core and one that has 
a radiative interior This is particularly interesting as the stel- 
lar parameters are almost the same for both stars, so their po- 
sitions in the CMD overlap within the usual observational un- 
certainties (see Table [TJ. The potential of this technique has 
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Fig. 7. Comparison for 1.3 Mq models using different mixing prescriptions: Ledoux (blue), Schwarzschild (red) and overshooting 
(green). Each row of models represents a value of the central hydrogen content (from top to bottom): Xc - 0.65, Xc - 0.45, 
Xc - 0.25, and Xc = 0.05. Left Panels: frequency ratios as defined in Eq. [14] Horizontal lines show a typical observable frequency 
range for each mixing prescription of radial orders n = 15-27. Central Panels: adiabatic sound speed showing the position of the 
density discontinuity. Right Panels: evolutionary tracks for the three selected prescriptions, and the position of the studied models 
is represented by a filled circle. Error bars have been plotted for reference in the models computed with the Schwarzschild criterion. 



only recently started to be exploited (iDeheuvels et al.l l2010t 
Ide Meulenaer et al.ll20ldl) . 

For more massive models, we show how this type of diagram 
can help us identifying different convective core sizes, evolution- 
ary stages, and therefore mixing prescriptions. Keeping in mind 
the results presented in Fig. |8] for a given mass value the dia- 
gram plotted in Fig.|9]can disentangle between different convec- 
tive core sizes, as tracks of a given mass and composition neatly 
separate according to the mixing prescription applied. With an 
accurate measurement of the large frequency separation, as ex- 
pected from the space missions, we can identify the size of the 
convective core and, as the mass and composition are known, the 
evolutionary stage. 

Nevertheless, mass values determined from a combination of 
asteroseismic measurements and modellin g are usually precise 
within 5-10% w hen metallicity is known ( Metcalfe et al.ll20Tol 
iGai et al.) 1201 Oh . In Table [1] we present the parameters of mod- 
els for each mixing prescription with a mass difference of up 
to ~8%. These are neatly distinguishable in Fig.|9]for the con- 
sidered error bars, showing that our method is able to break the 
degeneracy between mass value and mixing process. We stress 



again the fact that this analysis relies on known seismic and stel- 
lar parameters, in particular frequencies and metallicity. The ef- 
fective temperature differences for the mass range considered in 
Table [1] within a mixing prescription are ~150 K, a typical ob- 
servational error bar. Stars with a larger mass difference for this 
composition will have an effective temperature not compatible 
with the observations at the same large frequency separation, 
which allows us to exclude them from the analysis. 



As a closing remark, we notice that in Fig. |9]the track rep- 
resenting the evolution according to the Schwarzschild criterion 
shows abrupt jumps in its slope and mean values. This effect is 
a consequence of the convective process not properly taken into 
account at the edge of the core, and the abrupt variations corre- 
spond to the loops in the convective core profiles (see Sect. 13. lb . 
However, the same general behavior as in the other mixing pre- 
scription is present, so a mean value can be considered for com- 
parisons. 
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Fig. 8. Same as Fig.Qbut each row represents models with the same large frequency separation value, which is given in the right 
panel of each row (in fiHz). 



5. Discussion 

While investigating different convective zone boundary defini- 
tions, we have seen that using the Schwarzschild criterion for 
models with growing convective cores leads to inconsistencies 
in the exact location of these regions. The layers at the top of the 
core may suddenly be mixed into the central region increasing 
the central hydrogen content and developing loops in the convec- 
tive core boundary evolution (breathing pulses phenomenon). 
The question th at arises then is h o w to t reat these zones. In the 
study made by ICrowe & MitalasI (Il982h . the authors took into 
account these layers by mixing them until convective neutral- 
ity according to the Schwarzschild criterion was recovered. The 
convective core grows beyond its previously established bound- 
ary and a new supply of fresh hydrogen is added to the core, also 
increasing the timespan a star spends transforming hydrogen into 
helium. 

Another option for dealing with this breathing pulses is using 
convective overshoot as implemented in GARSTEC, producing 
a smooth convective core profile during the main sequence evo- 
lution. This prescription does not modify the thermal structure 
of the overshooting layers. When an adiabatic stratification of 
the convective zone is considered instead, small fluctuations at 
the edge of the core still appear but of smaller size than the ones 
presented in Fig.|2](M. Salaris, private communication). These 



fluctuations are not likely to affect the computation of adiabatic 
frequencies. 

Overshooting is an interesting topic by itself, and a detailed 
analysis of it goes beyond the scope of this paper However, we 
can mention two important aspects we think deserve attention. 
One is the temperature gradient in the overshooting region. We 
have assumed in our study that the temperature stratification in 
the radiative zone is not affected by the convective motions be- 
yond the formal boundary, while convective penetration assumes 
an adiabatic stratification in this layer Currently there is no con- 
sensus on the issue, and it is expected that numerical simulations 
and asteroseismology will cast some light on the topic (Godar^ 
2007). The other important aspect is the size of the overshoot- 
ing region. An exponential decay of the convective velocities is 
an alternative to an adiabatic extens ion to a fractio n of a scale 
height of the convective region (see IVenturalF2007l for details). 
However, numerical simulations are not yet abl e to reproduce 
the conditions of stellar interiors (lDintransll2009h . 

Regarding semiconvection, our calibrated value for the the 
overshooting efficiency washes out the molecular weight dis- 
continuity that produces semiconvective regions. Nevertheless, 
it has been suggested that the presence of a molecular weight 
barrier at the edg e of the core m ight reduce the amount of over- 
shooting from it (ICanutolll999h . It might be possible for a semi- 
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Table 1. Main properties of the models for a particular value of the large frequency separation. For each mass value, the structure of 
the computation for a given mixing prescription is given: Ledoux (Led), Schwarzschild (Sch) and overshooting (Ove). The columns 
are: effective temperature, gravity (logarithm), stellar age, central hydrogen content, size of the sound speed discontinuity at the 
edge of the convective core (Jump), position of the edge of the convective core in terms of acoustic radius (Aradcc), and the slope 
and mean values of the frequency ratios in the radial order n= 15-27. 
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Fig. 9. Slope versus mean value diagram of the frequency ra- 
tios in the assumed observable range throughout the main se- 
quence evolution. Top panel: 1 . 1 Mq models without overshoot- 
ing (solid red line) and with overshooting (dashed green line). 
Bottom panel: different mass values for the three considered 
mixing prescriptions: Ledoux (blue), Schwarzschild (red) and 
Overshooting (green). A representative 1-cr error bar is plotted 
on the top left. Models marked with a filled circle have the same 
large frequency separation value, and error bars correspond to 
3-cr uncertainties in the frequencies. See text for details. 



convective region to exist even with the presence of overshoot- 
ing. 

We have shown that the frequencies in the linear range are 
sensitive to the central regions of the star However, it is not 
straightforward to determine if the modes are reaching the con- 
vective core or not. The radial (/ - 0) modes propagate down to 
the center of the star, but the situation is different for the / - 1 
modes. Different estimations of the turning points can be made 
if the perturbation in the gravitational potential is neglected (the 
cowling app roximation. Cowling 1941) or taken into account 
dJeansI 19281) . In fact, it has been suggested that for stars with 
convective cores the dipolar modes reach the very center of the 



star (lTakatall2005Ll2006l) . This is an interesting topic where fur- 
ther development is surely needed. 

From the results presented in Sect. |4] we can say that it is 
possible to estimate the size of the homogeneously mixed zone 
using seismic data. In particular, we can discriminate between 
stars with and without a convective core. There are a number 
of assumptions behind this prediction. Besides the expected fre- 
quency data, we assume that there are determinations of stellar 
parameters from either spectroscopy or photometry, with an un- 
certainty in eff'ective temperature of ~150 K. Given these ob- 
servables we are only limited by the value of the stellar mass, 
which might be estimated within 5-10% precision for a given 
mixing presc ription u sing isochr ones and as t eroseismic scal ing 
relations ( Stel lo et al ]|2009; Basu etal 1 120101: iGai et d]l2010h . 

The frequency ratios seem to sufliciently isolate the signal 
arising from the stellar core allowing us to study the central con- 
ditions of a star. Another asteroseismi c tool aimed to study smal l 
convective cores was proposed by Cunha & Metcalf j (l2007l) . 
where a combination between small and large separations was 
constructed and analy zed using modes of angu lar degree / = 
0, 1 , 2, 3. In their work. lCunha & Metcalfel (12007 ) found that the 
sharp discontinuity formed at the edge of the convective core 
has a direct impact on the slope of their proposed quantity when 
plotted against frequency, showing that the absolute value of the 
slope increases as a function of the size of the chemical discon- 
tinuity (and therefore with decreasing amount of central hydro- 
gen content). Since observing I - 3 modes is very challeng iiig 
due to severe amplitude reduction. ICunha & Brandaol(l2010b ex- 
tended this study using the small frequency separations defined 
in Eqs. [12] and [13] and analy zing the changes in the slo pe, reach- 
ing the same conclusions as I Cunha & Metcaifg (120071) . 

The behavior with evolution is the same for the small fre- 
quency separation s, the frequency ratios, a nd the asteroseismic 
tool proposed by Cunha & Metcalf j (12007 ). Thus, our findings 
are applicable to all of them. We have shown in Sect. l4.3l that. al- 
though there is a relation between the size of the discontinuity in 
the adiabatic sound speed and the absolute slope of the ratios, the 
total extent of the convective core and the central hydrogen con- 
tent heavily influence the frequency ratios. Three models with 
the same large frequency separation and similar size of chemical 
discontinuity at the edge of the core have very different abso- 
lute slopes due to their different convective core sizes and cen- 
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tral hydrogen contents. On the other hand, models which barely 
change the size of the convective core during the hydrogen burn- 
ing phase show the mildest change in the absolute slope of the 
ratios. 

From the analysis performed in this paper, we have a clear 
theoretical understanding of the different contributions in the lin- 
ear range of the frequency domain and how to disentangle be- 
tween the effects of mass and type of mixing for realistic er- 
ror expectations in the observables. We have based our study 
of the frequency ratios assuming that the oscillations are pure 
p-modes, but this assumption breaks down when modes with 
mixe d p-mode and g-mode character appear (see lAerts et aP 
and references therein). These mixed modes carry infor- 
mation from the deep interior and are very sensitive to the core 
conditions, making them ideal probes for age caUbrations in 
sub-giant stars and mixing processes occurring during the main 
sequ ence (e.g 'Dziembowski & Pamvatnvkh"1991'; 'Straka et al.' 
12005.: .Deheuvels & Michel 2010- Metcalfe et al.. .2O10). For in- 
stance, they could be detected in main sequence F-type stars 
with shallow convective envelopes, where this avoided crossings 
could more easily reach the visible frequency range. However, 
different analysis tools from the ones presented in this paper 
must be developed to fully exploit their potential. 

6. Summary and conclusions 

We have presented results of our study of the influence of differ- 
ent convective boundary definitions and prescriptions of chemi- 
cal mixing on main sequence evolution of low-mass stars. This 
has important consequences for the overall evolution of low- 
mass stars. We have focused our study on the effects on stellar 
ages and CMD morphology of the existence and size of convec- 
tive cores. 

The age of a star at the end of the main sequence phase de- 
pends on the amount of fuel available for core burning. In the 
prescriptions we have considered, this is governed by the size of 
the homogeneously mixed zone in the center The implementa- 
tion of the Ledoux criterion for convection and no mixing be- 
yond the formal boundaries of the central convective region pro- 
duces the smallest size of convective core and consequently the 
shortest main sequence age for a given stellar mass and composi- 
tion. The time-dependent diffusive mixing in the semiconvective 
region increases the amount of hydrogen in the center, reach- 
ing the same convective core size as the Schwarzschild criterion 
for high enough mixing efficiencies. Convective overshoot can 
increase even further the extension of the mixed zone and the 
associated age value at the end of the main sequence. 

The appearance of a hook-like feature at the end of the cen- 
tral hydrogen burning phase can be shifted to higher mass values 
by the inclusion of the Ledoux criterion, as we have shown for a 
1.2 Mq model. This has relevant consequences when isochrones 
are used to reproduce the shape of the CMD of a stellar popu- 
lation given a chemical co mposition and a set of input physics 
dVandenBerg et al..,2007; Magic et al.ll2010 ). Previous methods 
of analysis do not allow to ensure the existence or absence of 
a semiconvective region in stellar interiors. Further development 
of numerical simulations and asteroseismic tools, such as the one 
prese nted in this paper, can be the key to solve this interesting 
issue (lBascou]|l2007l:lMigho et al.ll2008h . In the same Une, con- 
straining the existence and size of the central mixed region is 
another step towards a comprehensive understanding of mixing 
processes. 

We have focused our attention upon using frequency combi- 
nations sensitive to the convective core size and the evolutionary 



stage of the star to discriminate between the different mixing 
prescriptions. For this purpose, we have chosen the frequency 
ratios and analyzed their behavior throughout the main sequence 
evolution for different masses and mixing prescriptions. The ra- 
tios present a low-frequency domain which is governed by the 
global seismic properties of the star (the large frequency separa- 
tion). However, for higher frequencies the behavior of the ratios 
changes together with the inner conditions of the star. These are 
all high order low-degree modes and correspond to what we have 
termed the linear range. 

Using this linear range, our findings predict that we can 
discriminate between stars with and without convective cores, 
constrain the size of the homogeneously mixed central region 
and disentangle the effects from stellar mass. All of these pro- 
vided we use an adequate frequency combination to isolate the 
signal coming from the core. The method relies on our knowl- 
edge of the stellar composition, where large uncertainties can 
still be present. In order to fully assess this and other poten- 
tially important effects, such as helium abundance and micro- 
scopic diffusion, a hare and hounds exercise with extension to 
other metallicities and applications to real data is currently un- 
derway to further test the method (Silva Aguirre et al., in prepa- 
ration). Nevertheless, the extensive and highly accurate sets 
of o bservations obtained by the CoRoT and Kepler missions 
(e.g 'Michel et al.' '2008; 'Stello et al. " 20101: IChaplin et all 120101: 
IChristensen-Dalsgaard et al. 2010) open an exciting possibility 
for stellar physics: to apply differential analysis to stars with a 
common property, such as metallicity, and pin down the differ- 
ences arising from the inner layers of their structure. 
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